O QUE ESPERAR DO CALOR EM UMA PLACA ESTREITA ? 
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Resumo. Ilustramos com um simples exemplo de perturbagao de dommios, como a de- 
composigao espectral de um operador linear pode ser aplicada para determinar o com- 
portamento de uma familia de equagoes diferenciais parciais. 



1. iNTRODUgAO 

Na modelagem matematica de fenomenos naturais, e comum presumir que certas gran- 
dezas sao despreziVeis comparadas a outras. Por exemplo a profundidade de um oceano 
comparado a sua extensao. De todo modo, a grandeza desprezivel existe tanto quanto 
as outras, e a matematica leva isso em conta, atraves da nogao de limite: Se e existe no 
modelo, o que podemos dizer ao fazermos e — t- ? Em cada caso, explicar o significado que 
se da para "quando fazemos e — t- 0" , e explicar o que sera entendido por "desprezar e " . 

Muitos desses aspectos estao incorporados nos chamados problemas em dommios finos, 
sobre os quais [7j contempla um panorama substancial. 

O que consideramos neste trabalho exemplifica algumas ideias envolvidas no assunto, 
e, por isso mesmo, ap6ia-se em resultados importantes ja estabelecidos. De posse deles, 
o tratamento que damos ao exemplo vem com o objetivo de valorizar ferramentas basicas 
da Algebra Linear no tratamento de EDP's nao lineares. Elas, juntamente com o fato de 
podermos tornar preciso o significado de "fazer e — )• 0" , neste caso, permitirao darmos uma 
resposta quanto ao que esperamos ao fazer isso. 



2. MOTIVAgAO 

Imagine que queiramos estudar como o calor se propaga em placas retangulares de 
espessuras que vao tornando-se mais e mais finas. Fixando o comprimento das placas 
constante, digamos igual a 1, como prototipo tomamos o retangulo Q'^ := (0,1) x (0, e), 
e > 0, e consideramos e — )■ 0. 

Alem disso, para um modelo mais simplificado, supomos que as placas sao compostas de 
material homogeneo e estao termicamente isoladas. Isso significa que o fluxo de calor, sem 
fontes externas, e com distribuigao inicial de temperatura, Uq, supostamente conliecida em 
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cada ponto {x,y) do retangulo ^l^, e descrito pelas equagoes 

( K = ul^ + uly, em M+ x 

^ = 0, em M+ X 89." (Pe) 
^■^(0, •, •) = Uq, em ri*^, 

onde r]^ denota o campo vetorial (unitario) normal as laterals dil.^ da placa il.^ e u'' = 
u^{t,x,y) representa a temperatura no ponto {x,y) de fi*^ no Instante t > 0. Observe que 

— — it, x,y) = nos pontos de dQ'^, nos dlz que nao ha quaisquer perda ou ganho de calor 

atraves das laterals dQ^, em qualquer Instante t > 0. 
O operador de Laplace 

Au{x, y) = u^x{x, y) + Uyy{x, y), 

e um operador linear, e e na anallse de sens autovalores e autofungoes que seremos capazes 
de determlnar o comportamento llmlte do lluxo de calor em Q^, quando e — )■ 0. Note 
que, quando e — )■ 0, se colapsa no segmento (0,1), e assim suspeltamos que se tal 
comportamento llmlte exlstlr, devera ser regldo por um modelo unldlmenslonal. 

3. Analise Espectral 

3.1. Problema de autovalor. Para il*^ = (0,1) x (0, e) conslderemos o problema de 
autovalor 

—Au = Xu, em il.'' 

1^ = 0, em dn^. ^^'^ 

Para cada e > fixo, atraves do metodo de separagao de variavels (veja [4]), e facU obter 
que (todos) os val ores de A para os quals existe uma fungao nao nula, u = u{x,y), que 
satlsfaz a equagao (A^), sao 

'^mn = ^^("^^ + m,n = 0, 1, 2, • • • 

e que, alem dlsso, a respectlva autofungao assoclada a ^ e 

'^m,ni^^y) = cos{rmr x) cos{— y) . 

€ 

Munido da norma 

1 

2 



||n|| = yj \u{x,y)\ dxdy 

o espago das (classes de equivalencla das) fungoes cujo quadrado e Integravel a Lebesgue 
(veja por exemplo [8]), e complete e usualmente denotado por L'^{Q'^). Tambem, como esta 
norma e Induzlda pelo produto interno 

{u,v) = u{x,y)v{x,y)dxdy. 
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e um espago de Hilbert. Autofungoes, e u^^,^,, associadas a autovalores 

distintos, „ e X^, ^, respectivamente, sao ortogonais. De fato, integrando por partes, 
ve-se, que para quaisquer {m,n) ^ {m',n'), 



,n' TT 



(u^„,M^/„/)= / cos(m7r x) cos( y) cos(m' vr x) cos( y) dxdy 

, , , , , , / .niT n TT 
cos(m vr x) cos(m vr x) ax / cos( y)cos{ y)dy = u. 

Jo ^ ^ 



Alem disso, 



1 2 dxdy 



onde ao,o = 1, am,o = «o,n = e a^, 
Portanto as fungoes 



I cos(m7rx)|^(ix / | cos( y)\^dy 

Jo e 

: 2 para m, n = 1, 2, • • • . 



,n vr 



— cos(m7r x) cos( y), m,n = 0,1,2, 

^/e e 



constituem um sistema ortonormal em L^(Q'^). Como tambem e possivel mostrar, (veja 
[1]), que esse sistema e complete temos 



m,n=0 



(1) 



E facil ver para T : L^(r2'^) — ?• L'^iVL'') linear e continuo, que a convergencia das somas 

oo 

parciais em ([T|) implica em Tu = {u,w^^^)Tw^^^. 

m,n=0 

Entretanto, mesmo nao sendo —A um operador continuo de L^(n^), tambem temo^ 



J2 ^'^''^'rn,n) >'m,nW'm,n^ VMGDom(-A), 



(2) 



m,n=0 



onde Dom(-A) = {u e L'^{Vt') : ^ [{u,w'^^^) X'^^J^ < oo}. 

m,n=0 



3.2. Mudando de variaveis. Ao considerarmos a mudanga de variaveis 

: (xi,X2) ^ (x,y) := (xi,ex2) 



Ao leitor interessado sugerimos a segao sobre operadores compactos e auto-adjuntos em [T]. 
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que transf orm a a placa (fixa) Q = (0, 1) x (0, 1) na placa Q^, se u = u{x,y) e solugao do 
problema (A^), entao, pela regra da cadeia, v = u o ^ solugao do problema 



-Vxixi - \vx2X2 = , em O 
VxiVxi + -5^X2^X2 = 0, em 50, 



onde rj = {'i]xi,r]x2) denota o campo vetorial (unitario) normal as laterals da pla ca . 
Similarmente, se v = v{xi,X2) e solugao de d^eD , entao u = v o <I)^^ e solugao de [A^). 
Portanto a mudanga $e transforma solugoes do problema (A^) em solugoes do problema 
d^el ). Da mesma forma, as solugoes do problema (Pe) sao transformadas em solugoes do 
problema 



vl = <,xi + ^^X2X2' em M+ X 



^liVxi + ^<2^^2 = 0> em (90, 
v^{0, •, •) = Ug, em O, 



considerando-se = o 



Calculos diretos mostram que, para cada e > 0, os autovalores do operador 

1 



n 



2 



sao OS mesmos XL, „ = tt (m H — ^), m, n = 0, 1, 2, • • • , porem as respectivas autofungoes 
(normalizadas) agora independem de e e sao dadas por 

^{xi, X2) = am,n cos(m -K xi) cos(n7r X2), m, n = 0, 1, 2, • • • . 

Como antes, as fungoes fm,n) "^i^^ = 0, 1,2---, constituem um sistema ortonormal 
completo de L^(0), e 

00 

-Aen= ^ ('i^,W^,n>^m,n^^m,n> G Dom(-A,). 

m,n=0 

Ainda por separagao de variaveis, somos capazes de expressar a solugao v'' do problema 



( ) como 



v\t,xi,X2)= ^ e *^->" (T;^,'t;^^„) <,^„(3;i,X2). 



(3) 



m,n=0 



Escrevamos agora, para cada e > 0, os autovalores A,^ „ em ordem crescente (levando 
em conta suas multiplicidades 



< A^ < A^ < 
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Como, para cada natural k, existe um valor de parametro e^, que pode ser tornado de 
modo a que tt k < entao, sempre que < e < e^, 

Ai = 0, A2 = Xifi, A3 = A2^0' ■ ■ ■ ' ^k+l = ^fc,0- 

Portanto 

A^^vr2(n-1)2 := A°, 

para todo n = 1, 2, • • • . 

Da mesma forma, a sequencia de autofungoes v^, associada a sequencia A^, converge 
trivialmente (quando e — )■ 0), para a fungao 2:2) := \/2 cos(n7r xi), V n = 1, 2, • • • . 

Este comportamento dos autovalores do operador — A^, no qual o primeiro autovalor 
"associado a decomposigao na diregao X2" explode quando e — )• 0, e fundamental para a 
caracterizagao do limite para dP^D , pois ao reconhecermos e v^, n = 1, 2, • • • , observamos 
que sao, respectivamente, os autovalores e autofungoes (normalizadas) do seguinte problema 
unidimensional 

J -Vxx = Xv, em (0, 1) 

\ ^;,(0) = = 0. ^ 

Similarmente, o espago i^(0, 1) das fungoes definidas no intervalo (0, 1), com quadrado 
integravel a Lebesgue, e munido do produto interne 

{u,v) = / u{x)v{x)dx 
Jo 

e um espago de Hilbert, as fungoes constituem um sistema ortonormal completo, e 

£ ~ ^2 

k=l 

onde Dom(- — ) = {v G ^^(0, 1) : [{vA) >^k]' < 00}. 
^ k=l 

Tambem a solugao = v^{t, x) do problema unidimensional 



e dada na forma 



vt = Vxx, em M+ x (0, 1) 

Vx{t,{))=Vx{t,l)=^,t>Q (Po) 

v{S)r)=vl em (0,1). 

00 

{t,x) = Y,e-'>^l{vlvl)vl{x). (4) 

k=l 
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4. REDUgAO DA DIMENSAO 

Dada uma fungao u : (0, 1) — t- M, consideramos n : J7 — t- M definida por u{x, y) = 
u{x). Assim podemos ver o espaQO -L^(0, 1) como um subespago de L^(r2), e considerar o 

operador — atuando neste subespago. Na sequencia usaremos sem qualquer mengao 
esta identificagao. 

Com a notagao acima, os vetores constituem um sistema ortonormal completo de 
L2(0,1), e 

oo 
k=l 

Portanto, da convergencia espectral obtida na Segao El se a sequencia de dados iniciais 
Vq G converge para Vq em L^(il), entao 

oo oo 
k=l k=l 

uniformemente para t em intervalos compactos de Sobre estes fatos indicamos ao leitor 
as monografias [3] e [5]. 

Portanto, apos essa nossa analise, se perguntados fossemos : 

- O que esperar do calor em uma placa estreita ? 
de pronto responderiamos: 

- Esperamos que ele se remeta ao caso unidimensional. 
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